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ABSTRACT 

We build and test Parker-wind models to apply to observations of large-scale (~ lOOpc) outflows from Active 
Galactic Nuclei (AGNs). These models include detailed photoionization simulations, the observed radially 
varying mass profile, adiabatic cooling, and approximations for clouds dragged along in the wind and the 
interaction of the wind with the circumnuclear ISM of the galaxy. We test this model against recent HST/STIS 
observations of [O III] emission-line kinematics (in particular, we test against those observed in NGC 4151, 
but approximately the same kinematics is observed in NGC 1068 and Mrk 3) to constrain the viability of large- 
scale thermal winds in AGNs. We find that adiabatic cooling dominates in these outflows, decelerating Parker 
winds on large scales, making them highly unlikely as explanations of the observed kinematics. 
Subject headings: galaxies: active — galaxies: Seyfert — galaxies: emission lines — hydrodynamics — 
radiative transfer — galaxies: individual (NGC 4151, NGC 1068, Mrk 3) 



1. INTRODUCTION 

Outflows are ubiquitous in astrophysics: they rather para- 
doxically accompany collapse processes, whether the collapse 
of a molecular cloud core to form a young stellar object, the 
collapse of matter onto a galaxy's central black hole, or the 
collapse of stars into white dwarfs, neutron stars, or black 
holes. Outflows are important to study for their links to those 
basic collapse and accretion processes, but also because un- 
derstanding outflows is crucial to studying the impact and 
feedback of those winds on the surrounding medium. 

Outflows from AGNs are interesting for both of the reasons 
given above. While there is widespread agreement that 
AGNs function via accretion from subpar sec-scale disks 
onto super-m assive black holes (see, e.g., Peterson 1997; 
iRrolikl fl 999ft . the detailed physical processes controlling, 
for instance, the generation of the continuum, obscuration 
of AGNs, transferring matter from large (~ 1 kpc) to small 
scales, as well as the launching mechanisms for outflows, 
are still not well understood. Most importantly here, while 
^50% of AGNs show evidence for outflowing U V and 
X-ray absorption jRevnolds et alJll997t Ideorge et all Fl 998: 
ICrenshaw et all Il999t iKrissI 120011) . the physical mecha- 
nisms powering those outflows remain unclear. Magnetic 
fields are a leading candidate for launc hing winds in both 
AGNs and young stellar objects ( Blandford & Pavne 
1 9821 lEmmering. Blandford & Shlosmanl [ 1992]: 

Contopolous & Loyelacd 119941: iBottorff etail 1 19971: 
Bottorff & Ferlanl 120001: IBottorff. Kori sta & Sh losmanl 
2000). In addition, given the high-flux radiation field in 
some AGNs (especially the high-luminosity quasars), it is 
also natura l and useful to examine radiative acceleration of 
outflows JMcKee & Tarterl 11 9751: lShieldsl[T977t ITckd 119771: 
Shlosman. Vitello & Shayivl Il985t lArav. Li & Begelman 
19941: iMurrav et alJ Il995t IChiang & Murravl |1996[ 
Proga. Stone & Kallmanl 12000 IChelouche & Netzeil boOll 
2003a b; Proga & Kallman 2004). Some groups have 
also combined these two mechanisms, addressing the 
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import ance of radiative ac c eleration within ma gnetic 
winds JKonigl & Kartiel 119941 Jde Kopl & Begelmanl 1 19951: 
lEverett. Konigl. & Aravll2002l: lEverettl 120051) . Apart from 
these two models, thermal wind models have also been 
proposed; several papers have exam i ned X-ray heated winds 
( Bege lman. McKee & Shields! 1 19831: IKrolik & Vrtileld Il984t 
Balsar a & Krolikl [19931 iKrolik & Krissl 119951 IWoods et alJ 
1 19961 lKrohklll997l) . IChelouche & NetzeJ fc005l) examined 
the possibility of thermally-driven (Parker) winds and the 
ability of those winds to explain observations of low velocity 
(v ~ 500 km s" 1 ) X-ray absorption features. In this paper, we 
also examine the possibility of Parker winds, but in a different 
context; we test these wind models against spatially resolved 
observations of the kinematics of forbidden-emission lines 
(here, [O III]) in local, low-luminosity AGNs. 

Forbidden-line emission is very useful for testing models of 
AGN outflows and their interaction with surrounding gas be- 
cause such emission is seen at relatively large distances from 
the nucleus (of order 100 pc). A large distance for the emitting 
gas is consistent with the observed narrow (low velocity dis- 
persion) emission lines; this region is therefore known as the 
Narrow Line Region (NLR). The large distance to the NLR is 
important because whereas many of the physical processes in 
AGNs occur on scales much too small to directly observe and 
resolve, the NLR can be and is resolved in some AGNs. 

Resolving the NLR is largely a product of the high reso- 
lution observations of the Hubble Space Telescope, although 
early ground-based studies of NGC 4151 revealed t he pos- 
sibilit y of large scale accelera tion a s early as 19 90 (Schulz 
H990t TMoore. C ohen & MarcvlT996L iMoore & Cohen 1996). 
More recent spectra taken with the HST's Space Tele- 
scope Imaging Spectrograph (STIS), however, allow high 
spatial and spectral resolution studies of the NLR, offer- 
ing data sufficient to constrain photoionization and dy- 
nam ical mode l s of the outskirts of AGNs. In particu- 
lar, iDas et all J2005I) have mapped the kinematics of up 
to three different flux components of [O III] emission in 
NGC 4151 in five different slits across the observed NLR. 
In addition, previous studies of the NLR of NGC 4151 
( Wing e etalJll999t iNelson et alJl2000t iKraemer et al]l200fl> 
and of the NLR of other AGN (e .g., Crenshaw & Kraemer 
120051 IKraemer & Crenshaw! l2000l iRuiz et alJ 12005ft have 
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found similar kinematics. The most striking result from 
these studies is the apparent near-linear increases in veloc- 
ity (v oc r) over scales from 10 pc to 100 pc and then an 
apparent linear decrease in velocity at larger scales; fur- 
ther, this trend has been obser ved on similar sc ales in three 
different objects: NGC 415 1 <Das et al]|2005l). NGC 10 68 
dCrenshaw & Kraemerl2000l) . and Mrk 3 JRuiz etalJ2001l) . 

These observations may help constrain the driving forces 
behind large-scale AGN winds. At first glance, both radiative 
acceleration and magnetic acceleration are problematic with 
respect to these observations: both mechanisms accelerate gas 
to its escape velocity on length scales of order the launching 
radius. If the outflow is launche d from an accreti on disk, the 
launching scale must be < 1 pc (Goodman 2003), and there- 
fore both radiative and magnetic winds would be expected to 
reach their terminal velocities by approximately 10 pc from 
the central source. Another possibility for the driving is ther- 
mal, or Parker, winds: could such thermal winds accelerate 
gas slowly enough to explain the observations? This seems 
plausible, especially if we allow the [O III] emission to orig- 
inate in clouds dragged along with the wind, which would 
further slow the observed acceleration. We can then ask 
what constraints these observations of v oc r outflows place 
on thermal wind models. Also, just as important: in the three 
cases mentioned above, what decelerates the gas? Finally, 
such models ma y help constrain the geometry of the NLR 
dDas et alJ2005h . 

To address these questions we must build detailed kine- 
matic and photoionization models of the core of AGNs and 
their winds. In this paper, we build such a thermal wind 
model (in ^2), adapted to considerations of large-scale flows 
in AGNs, and then apply that model to HST/STIS observa- 
tions of NGC 4151. The particular spectral energy distribu- 
tion and M(r) profile in NGC 4151 are introduced in |3] The 
particular mass profile that we adopt allows for a good first- 
order approximation to the wind as roughly isothermal, so 
isothermal wind models are applied to NLR observations of 
NGC 415 1 in ^] Isothermality in AGN winds is then tested in 
the particular case of these thermal winds (and tested against 
variations in the input spectrum) in ^5] We find that thermal 
wind models face very stringent constraints that make it very 
unlikely that they can explain the observed NLR kinematics, 
due chiefly to the effects of adiabatic cooling which ensures 
that the winds are not isothermal, and in fact cool on large 
scales (r > 5 pc). 

2. PARKER WINDS 

We first review the basic Parker (solar) wind model f ^2.H . 
and then introduce changes to that model that must be in- 
cluded when hypothesizing a large-scale Parker wind at the 
center of an active galaxy: the increase in enclosed mass as a 
function of radius, adiabatic cooling for outflows at such large 
distances, and finally the effects of wind/cloud and wind/ISM 
drag. 

2.1. Parker Winds Applied to AGNs 

A Parker wind can be thought of as an extended ther- 
mal wind. To explain this terminology, first consider the 
simplest picture of a thermal wind, where gas at its launch 
point is heated to high temperature such that that its sound 
speed is greater than the escape speed (c s > v esC ape)- This 
requirement can be recast to ask if the gas temperature 
is greater than the "escape temperature" r e = GMfi/Rk 
(Begel maiTMcKee & Shield! 198 3). where M is the mass of 



the central object, /i is the mean mass per particle, R is the 
launching radius, and k is Boltzmann's constant. An inher- 
ent difficulty with this picture, in relation to the NLR data 
discussed in the Introduction, is that there is no extended ac- 
celeration as is apparently observed. 

A Parker wind also depends on thermal effects, but with 
the major difference that the gas is not heated at the base 
of the wind to have c s > v e scape- The wind is therefore more 
like an extended atmosphere, where, at the base of the wind, 
Cs *C v escape . The key difference, though, between a Parker 
wind and an atmosphere is that if such strongly bound gas is 
heated as a function of height such that the temperature drops 
more slowly than M(r) / r, where M(r) is the enclosed mass 
and r is the spherical radius, the atmosphere must expand su- 
personically (explained in more detail later in this section; see 
Parker! Il965l) . This leads to rather slow acceleration, com- 
pared, for instance, to radiatively-driven and magnetically- 
driven winds. 

As mentioned previously, therefore, this model is there- 
fore appealing for the NLR of AGNs. This is espe- 
cially true of NGC 4151, as photoionization analyses have 
shown that the c entral AGN continuum dominates the heat- 
ing of the NLR ^Alexander et al.lfl99^ iNelson et all 1200(1 
iKraemer et ail 120001 12001 and references therein). These 
photoionization analyses also show that shock heating (which 
could, for example, result from an interaction w ith a colli- 
mating medium, e.g. Brighenti & Mathews 2006) is also not 
likely to be important. (We note, though, that Kraemer et 
al. 2001 find the possibility of shock heating near the start 
of deceleration of the wind at approximately r ~ 100 pc.) 
It is also possible in some AGNs for the central radio jet 
to affect the NLR by increasing the NLR line widths (e.g., 
[ Nelson & Whittlell996tlC aDetti. Macchetto & Lattanzil 19971: 
Steffe n et al.ll997t |Axon et al. 1998) but that does not appear 
to be the case for NGC 4151 JDas et alJ2005l) . 

To examine this possibility in more detail, we start with the 
basic equation for an isothermal Parker wind: 



dv 
dr 



2c 



GM 



•0-5) 



(1) 



where v is the one-dimensional radial velocity, and r is the 
spherical radius. 

Two important insights can be seen from examining Equa- 
tion n closely. First, we note the importance of the "critical 
point" in this equation. As stated earlier, gas at the base of a 
Parker wind is gravitationally bound, so that GM/r 2 ^> c 2 Jr. 
Therefore, the numerator of Equation ^ must be negative at 
the base of the wind; since the denominator is also nega- 
tive (v < c s ), a positive acceleration results. However, as the 
wind moves outward and v increases, eventually v = c s and 
the denominator of Equation ^ becomes zero. In order for 
the Parker wind equation to be physically meaningful, the 
numerator must also be equal to zero at the same time the 
denominator becomes zero, so 2c 2 /r must equal GM/r 2 (or 
r = Critical = GM/2c 2 ) must be true when v = c s . This require- 
ment sets the location of the critical point: this is the (r,v) pair 
that the solution must "thread" so that the wind solution can 
describe a flow that continuously accelerates. In fact, the solu- 
tion that passes through this critical point is the only solution 
that allows an outflow to large distances and to speeds above 
the escape velocity. Physically, this critical point represents 
the last point downstream in the outflow where the boundary 
conditions at the disk can communicate with and therefore af- 
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feet the wind. An example of such a solution is shown by 
the dashed line in Figure [l] which threads the critical point, 
indicated by the circle. Other solutions (such as the "breeze" 
solutions that have dv/dr = below the position of the criti- 
cal point in Fig.^ achieving smaller peak velocities than the 
critical point velocity and then decelerating) that do not result 
in a viable outflow are also shown. 




FIG. 1 . — A variety of solutions to the Parker wind equation (Eq.0 for an 
isothermal wind with T = 5 X 10 6 K, M. = 1.3 X 10 7 Mq, and Haunch = 0.1 pc. 
The dashed line shows the solution that passes from low velocity at the launch 
point, through the critical point, to high velocity at large distance. This solu- 
tion, passing through the critical point, is the only steady-state solution that 
yields a viable wind model. Other solutions, not passing through the critical 
point, are shown for comparison. For completeness, we also show that it is 
possible to have an inflow solution that passes through the critical point: gas 
at high velocity (v 2> c s ) at large distance (r ^> r cru i C ai) that slows as it falls 
to the central object; this is known as Bondi-Hoyle accretion I Bondi & Hovld 
1944). 



Second, in order to pass through that critical point, there is a 
further constraint on the temperature structure of the outflow, 
which was mentioned briefly above. In order for the solution 
to pass through the critical point and continuously accelerate, 
the temperature profile, T(r) must decrease with radius less 
quickly than M(r) jr. To show this, we recast the Parker wind 
equation into a form that shows directly the comparison be- 
tween the pressure gradient and gravity: 



dv 

r — 
dr 



dr 



GM 



(2) 



This equation is simply a more general form (not assuming 
isothermality) of Equation ^ which will allow us to investi- 
gate non-isothermal wind models. The same conditions that 
applied to the numerator of Equation[2apply to the right-hand 
side of Equation|2] namely, it will also have to transition from 
negative to positive values at the critical point. 

Now, instead of assuming an isothermal wind (as in Eq.Q, 
here we set T = T(r) = 7o(r/ro) Q , but retain the assumption of 
a central point mass, M. After substituting this into Equation|3 
and simplifying, we find: 



dr 



— -(a- 

r a v 



-2)/*- 1 =- 



GM 

T 



(3) 



where c s o is the sound speed at the base of the wind. One can 
see that in order for the pressure term to dominate gravity, a 



must be greater than -1. To illustrate this, consider the case of 
a = -1 in Equation[3] 



dr 



3c lo r o 



GM 



(4) 



At the base of the wind, the right-hand side of the equation 
will read 3cJ /ro — GM/r^. Since, at the base of a Parker 

wind, Cj <C GM/ro, this term will be negative at the base 
of the wind, and since both the pressure gradient term (the 
first term on the right-hand side) and the gravitational term 
decrease as r 2 , the right-hand side of the equation will never 
be positive. In that case, the wind will not accelerate to escape 
velocity. Therefore, in order for the wind to accelerate in this 
case (where M is constant with radius), a must be greater then 
-1. One can derive the same result anoth er way : if a < — 1, the 
result is a static atmosphere (see Parker 1963j). 

For large-scale winds (r > 1 pc) in AGNs, though, M{r) 
increases with radius. Including the possibility of a spatially- 
varying M{r) in our requirement on the temperature profile 
T(r), above, r a must decrease less quickly than M(r) /r in or- 
der to accelerate the wind to large distances. Assuming an 
isothermal sphere to model the gravitational potential of the 
circumnuclear matter, we will find that (for NGC 4151) the 
black hole stops dominating the gravitational potential at ap- 
proximately 1 pc. Beyond that distance, M(r) oc r. In this 
regime, assuming a simple M(r) = Mo(r/ro) for now, Equa- 
tion[3]becomes: 



dv 
dr 



L .v.O 



(a-2)r° 



GM 
ror 



(5) 



Equation[5]then shows that we must have a > in order to get 
large-scale acceleration in a thermal wind. In other words, the 
wind must have T(r) at least slightly increasing with radius in 
order for a Parker wind to accelerate. 

In order to check this constraint on T(r), we must ensure 
that we address all of the heat ing a nd cooling components in 
an AGN wind accurately. In 33.11 we will use Cloudy (Ver- 
sion 05.07.06; Ferla nd et alJl!998f> to simulate the photoion- 
ization heating, although Cloudy by itself cannot include the 
effects of adiabatic cooling in the context of our thermal wind 
and wind geometry (it can calculate adiabatic cooling in the 
context of its own internal wind models, however). This im- 
portant cooling source will be considered next. 

2.2. Adiabatic Cooling 

We will run photoionization simulations at regular 
positions along the flow, but those separate photoion- 
ization simulations will not have any knowledge about 
the kinematics of the outflow. This would problem- 
atic for these large-scale flows where adiabatic cooling 
is important (see, e.g. IRu den. Glassgold & Shu 1990; 
[ Frank. Norieg^ Crespo & Bali ckl 1 19921 ISafieiT 119931: 
iShangetalJ I2002L IChelouche & Netzerl 120051) . so we 
calculate the magnitude of the adiabatic cooling outside of 
Cloudy and include it as input to Cloudy. The adiabatic 
cooling is simply calculated from: 



dq de d ( 1 
— = — +P— - 

dt dt dt \p 



(6) 



where q is the heat per unit mass, e is the internal energy per 
unit mass, and P d(l/p) is the work per unit mass done by 
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the expanding gas. After some algebra, the change in internal 
energy is given by: 



d (3P\ _Pdp 
^ dt \2 p ) p dt 



(7) 



Finally, using M out = 4irr 2 pv to eliminate p from the above 
expression, we find that the adiabatic cooling is given by: 



w 2v a \ i 

adiabatic = pC s [ — + - ) ergS COl S 



(8) 



where p is the mass density of the outflow, c s is the speed of 
sound, v is the radial velocity, r is the distance from the central 
source, and a is the acceleration in the wind. This particular 
expression and the above-noted units are those necessary to 
supply the correct cooling rate to Cloudy ; as such, the units on 
this ex pression are different from, e.g., Dal garno & McCravl 
(1972). 

For all of the wind solutions found in this work, a/ v is at 
least three orders of magnitude smaller than 2v/r, so the a/v 
term is safely neglected. 

2.3. Clouds Dragged Via Ram Pressure 

We will model the [O III] emission as coming from clouds 
injected into the wind and accelerated by the ram pressure of 
the wind on those clouds. The drag forc e for each individual 
cloud is computed via (e.g., Smi th! 19841) : 



rag.cloud — Pwind(^wind ^cloud) -^cloud 



(9) 



where p w i n d is the mass density of the wind, v w i„d and v c i ou d 
denotes the velocity of the wind and clouds, and A c i ou d repre- 
sents the cross sectional area of a cloud. Assuming a spherical 
cloud with mass M c i ou d = (4/3)7rr 3 p c i ou d, the acceleration for 
each cloud is 



^drag, cloud " 



Pwind l(v w ind- Acioud) 



Pcloud 



^wind 3 , \2 

a7 T v^wind ^cloud/ 

^VH, cloud 4 

(10) 

where R c i on d is the radius of an individual cloud. The ra- 
tio «wind/^VH,cioud is a free parameter in our models. We 
use this parameter for all clouds, although it would be more 
physically plausible to imagine a range of densities for the 
clouds intercepting the large-scale thermal wind. This would 
lead to a range of velocities, which might explain the large 
(~ 500 km s" 1 ) [O III] line widths observ ed in this objec t 
jDas et alJl2005l) and seen in many AGNs llRice et all 2006). 
The model presented here solely addresses trends in the ob- 
served centroid of the lines, however, and so only a single 
density ratio is used. 

We assume that the drag of the clouds on the wind is negli- 
gible compared to the acceleration of the Parker wind. 

2.4. Decelerating the Parker wind via ISM drag 

The wind is decelerated in much the same way that the 
clouds are accelerated: ram pressure from the impact on an 
external medium, which we assume is also broken up into 
clouds, and is also assumed stationary. Much like before, 
then, the drag force is given by: 



fl drag,ISM clouds — 



Pvj'mA 



?-V 2 
4 wind 



^wind 



PlSM cloud ^ISM i 



A^H.ISM ( 



(ID 

As the wind decelerates, the embedded clouds decelerate 
as well, as v w ; n d becomes smaller than v c i ou d in the cloud drag 



equation, Eq[K)| Again, nwind/A^H.iSM cloud is a free parameters 
in the models we run. 

The wind decelerates under this drag force until the wind 
becomes subsonic, at which point the integration of the equa- 
tions of motion ceases. However, v < c„ at such large dis- 
tances that it does not affect our ability to fit the observed 
kinematics. 

3. APPLYING THE PARKER WIND MODEL TO NGC 4151 

The discussion in ^outlined the basic physics of the Parker 
wind. Next, the inputs to that model must be defined from the 
observations of, in this case, NGC 4151. Most importantly for 
the consideration of heating, the centra l spe ctral energy distri- 
bution must be realistically modeled ( 33. 1> an d the particular 
M(r) profile for NGC 4151 must be included ( CHI . 

3.1. Modeling the Central Spectral Energy Distribution 

An accurate central spectral energy distribution (SED) is 
very important to understanding the heating and subsequent 
acceleration of Parker winds in NGC 4151. To model the 
S ED, we use the central continuum of this object as modeled 
by Kraeme r et alJ (12000): this continuum is then used as input 
for the Cloudy photoionization simulations. The continuum 
has two components: an intrinsic SED and absorption from 
surrounding gas. The specification for the intrinsic, unab- 
solved SED is given in Tabled the first and last values of a in 
the table are relatively unimportant, and are specified simply 
to cut the spectrum off at low and high energie s, whil e the in- 
termediate values of a are from Kraemer et al. (2000). To this 
simple power law continuum, we add absorption via an X-ray 
warm absorb er and UV absorber ( components "X-High" and 
"DEa" from Kraemer et al. 2000). The "X-High" absorber 
has a total hydrogen column density of A^h = 10 22 5 cm" 2 and 
an ionization parameter of U = «ion/«H = Li 0n /A / Kr 2 nnc = 1 .0 
(where L lon is the ionizing luminosity for E > 13.6 eV), 
while the "DEa" absorber has a total hydrogen column of 
A^h = 10 19 75 cm" 2 with U = 10~ 3 (using the c olumn for the 
"northeast" absorber from Kraemer et al. 2000); the metallic - 
ities listed in Kraemer et al. (2000) are used for these absorp- 
tion systems and therefore used for all of the photoionization 
simulations presented here f or self -consistency. (The metal- 
licities from iKraemer et all (|2P00) includes only those met- 
als with abundances greater than ~ 10~ 5 of hydrogen, and for 
those metals, differs only by 0.2 dex from Cloudy's default 
solar metallicity. Using Cloudy's default solar metallicity for 
the remaining calculations result s in o nly insignificant differ- 
ences in the T(r) calculations in 35.11 ) The absolute luminos- 
ity was set to 10 43 56 erg s" 1 in the range of 13 .6 eV to 1 keV 
to mat ch the spectrum displayed in Figure 2 of lKraemer et alJ 
(2000). The resultant continuum is then compared to a range 
of observations in Figure 13 The generally good agreement 
between the model continuum and the observations is essen- 
tial for accurately modeling thermal winds in this source. 

It is important to note, of course, that the luminosity of 
NGC 4151 has been observed to vary by factors of a few 
over the past decade (see, e.g., Kraem er" et alJl2005l where 
the continuum level was seen to be roughly 2.5 times brighter 
in the 2-10 keV band in 2002 vs. observations in 2000). The 
continuum that we have assumed here fi ts reasonably well to 
the observations of Edelson et a DU996), taken in December 
1993 when NGC 4151 was near its peak historical brightness, 
so tests with this continuum should yield good estimates for 
maximum heating of thermal winds. We will also ch eck t he 
dependence of the results on the central continuum in 35.31 
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TABLE 1 

Parameters for NGC 415 l's intrinsic, unabsorbed 
SED 



in £end Power Law ol\ F v oc is a 

1.36 x l(r 7 eV 1.24 x 1CT 3 eV 2.5 

1.24 xl(r 3 eV 13.6 eV -1.0 

13.6 eV IkeV -1.4 

1 keV 100 keV -0.5 

100 keV 100 MeV -2.86 



: Transmitted Contiuum 

+ Edelson et al. (1996) 
r x Kriss et al. (1992, 1995) 
: o NED 




Rest Frequency [Hz] 



FIG. 2. — Comparison of the assumed central SED for NGC 4151 (defined 
in TablefTlwith added intrinsic absorption) with observations. The solid line 
shows the SED, modeled by Cloudy, after absorption by components "X- 
High" and "DEa" from Kraemer et al. 1 2000). The red data points show the 
optical, UV, X-ray and 7-ray results from (Edelson et al. 1 1996), the green 
data points show the rang e in UV brightness observed during HUT pointings 
iKriss et al. 1992, 1996), and the blue data points show the wide array of 
spectral data in NASA's Extragalactic Database. 



We next calculate the Compton temperature for this gas, 
using the above SED and luminosity to calculate the gas 
temperature at a range of values for the gas ionization pa- 
rameter, U . The resultant plot of equilibrium gas tempera- 
tures (as in Rrolik, McKee & Tarter 1981) is shown Figure^ 
which also shows that T e .compton = 3.3 x 10 7 K. Note, however, 
that at extremely high ionization parameters and temperatures 
(U/T e >Q.l;U> 2.5 x 10 6 and T e > 3.3 x 10 7 K), the equilib- 
rium temperature starts increasing again. This is due to stim- 
ulated Compton heating of the gas, which occurs only at very 
high ionization parameters: when the radiation field is very 
intense compared to the gas density (when the electrons are 
only an "impurity" in the flux of photons), the stochastic elec- 
tric field of the high-intensity photon field modifies the elec- 
trons' response to the photons, yielding higher h eating rates 
JLevich & Sunvaevll97(l Blandfor dTl973llWilsonll982i) . We 
include this regime in the equilibrium curves shown here for 
completeness, but note that even without adiabatic cooling, 
the ionization parameters in the Parker winds explored here 
never reach the ionization levels where this extreme heating 
dominates. In addition, including adiabatic cooling results in 
significant decreases i n this maximum temperature (as already 
noted, for instance, bv lChelouche & NetzeJ2005l) . Given the 
SED and equilibrium curve, we can ask if the simplest thermal 
winds could explain the NLR velocities in NGC 4151. For a 
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FIG. 3. — The equilibrium "S" curve for the central absorbed SED in Fig- 
urel2l This plot shows that the Compton temperature of gas near the central 
source in NGC 4151 is 3.3 X 10 7 K. Note the upturn in T e at U/T e > 0.1 due 
to stimulated Compton scattering. 



range of initial wind densities, we show the thermal speed 
and escape speeds for winds launched from 0.1 and 1 par- 
sec in Figure|4] For winds launched at 0. 1 pc, the gas is not 
heated to sufficient temperatures to escape the gravitational 
potential. Winds launched at r > 1 pc have sound speeds high 
enough that gas can be evaporated off the disk and escape 
the gravitational potential to approximately the correct veloc- 
ities, although such a simple evaporation model would not 
explain the large-scale acceleration seen in NGC 4151. How- 
ever, Figure |4] does show that high density winds will have 
c s o < C Ves cape, so that Parker winds may be viable solutions 

(see am . 




FIG. 4. — The sound speed, c s , vs »ih, the number density, at the disk, 
using the SED shown in Fig. [2] Only densities for which c s > v esc ape can 
gas be "evaporated" off of the disk. This figure shows that such evaporation 
only happens for large distances (r > 1 pc) and very small densities («h ^ 
3 X 10 3 cm -3 ). This gas is not appreciably accelerated at large distances as is 
observed, so this model cannot explain the observed velocities in NGC 4151, 
but it does hint that more detailed thermal winds might be viable models, 
since approximately the correct velocities can be achieved. 



3.2. Radially Increasing M(r) 
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As mentioned in 32.11 the increase in enclosed mass with 
radius is crucial to modeling thermal winds in AGNs (and 
has been c onsidered without photoionization simulations by 
iKrolik & Vrtilekl[l984l) . Observationally, the mass distribu- 
tion within NGC 4151 is not well understood: the only data 
that exists are estimates of the black hole mass, and estimates 
of the galaxy's potential at approximately 700 pc. As a first 
order approximation, we assume an isothermal sphere mass 
profile and so therefore connect the inferred masses with the 
following functional form: 



M(r) = M. + '■ 



Ira 1 



(12) 



where a(r) is the observed velocity dispersion at radius r and 
M(r) is the enclosed mass within radius r. 

The c entral black hol e mass is M. = (1.33 ± 0.46) x 
10 7 Mq dPetersonll20 04). This mass must dominate at small 
radii (with this M(r) profile, at r < 2 pc). At large radii, the 
enclosed mass can be derived from HI 21cm line studies; we 
find the best current constraint to be a(r = 10" = 690pc) = 
122 km s" 1 (Mundel l et aljfl99^) . These two data points de- 
fine our enclosed m ass f unction, M(r). 

As discussed in 32.11 M(r) oc r at large scales in AGN is 
very important for Parker wind models, as it sets stringent 
constraints on the T(r) profile required for the Parker wind to 
be launched to infinity. This M(r) profile also changes the po- 
sition of the critical point in the wind. Inserting Eauationll2l 
for M(r) into the Parker wind equation (Eq0 and solving for 
the critical point condition, we find 



'critical — 



GMn 



2(cj-a 2 ) 



(13) 



In the limit where Mbh is dominant (a = 0), this equation re- 
duces as expected to the standard expression for the critical 
point, r cl ; t i ca i = GM Bli /2c 2 . Further, the acceleration at the 
critical point is also modified from the standard Parker-wind 
expression because of the inclusion of the M(r) function, and 
is given by 



dvY 

dr )r=r a 



(c 2 s -a 2 f 
(GM B11 ) 2 



(14) 



These equations allow for integration of the wind equations 
starting from the critical point. 

4. ISOTHERMAL PARKER WIND MODELS 

To test the Parker wind model against observations, we first 
numerically solve the isothermal Parker wind model, integrat- 
ing from the wind's critical point to the disk and then from the 
critical point out to where the wind decelerates to v < c s . That 
ID solution to the Parker wind is then mapped to a biconical 
outflow (with geometrical parameters largely similar to those 
in Das et al.) and then compared to the kinematic data. Our 
geometry is depicted in Figure |5j the "near" component of 
the bicone, for NGC 4151, is close to being in the plane of the 
sky, while the "far" component of the bicone is close to being 
out of the plane of the sky. We will use this terminology in 
refering to the different kinematic components of the bicone 
that are used to explain the spread of observed velocities in 
our dynamical models. 

Given these assumptions about the outflow geometry, the 
best-fit solution to the data (only for th e "north-east" outflow 
in Slit 1 of the data in Das et al. 20051 as we adopt the geom- 
etry from Das et al. and do not fit data from other slits) for 




'near" 
bicone 

(nearly in the plar. 



observer 



(nearly perpendicular to the plar 



Wind Launching Point 



FIG. 5. — A diagram of the biconical geometry for the observed NLR 
outflow in NGC 4151. The observer, at right, sees the outflow tilted at an 
inclination angle of 45° to the line of sight. This yields to halves to the 
bicone with somewhat distinct kinematics: a "near" bicone that is nearly in 
the plane of the sky, and a "far" bicone with is closer to being perpendicular 
to the plane of the sky, or nearly along the line of sight of the observer. 



the case of clouds being dragged in the wind is shown in Fig- 
ure [6] where the velocities displayed are the velocities of the 
embedded clouds as a function of radius along the biconical 
outflow. This model does a reasonably good job of encom- 
passing the observed velocity data points. For completeness, 
the thermal wind velocities (the velocities in the continuous 
wind itself and not the embedded clouds shown in Fig.^Ji are 
displayed in Figure [7] The parameters for this fit as well as 
"error bars" are given in Table|2] The most important param- 
eter from this table and from the isothermal wind models is 
the required temperature, T e = 3 x 10 6 K. 

Note that the error bars in Table [2] are only very approxi- 
mate quantities that result from manually varying each param- 
eter individually, and inspecting a range of fits to the data and 
picking those fits and ranges of parameter values that seem to 
best represent the data. We are aware of no method to cal- 
culate a % 2 -like figure-of-merit to quantitatively analyze such 
data, but present these approximate error bars to communicate 
the sensitivity of the fit to various parameters. Most impor- 
tantly, these error bars show that the isothermal wind temper- 
ature is fairly well constrained, given the other parameters in 
the fit. Also, note that the initial radius is very much uncon- 
strained; the only constraint we have here is that ro < 0.8 pc, 
which comes from the requirement that ro < r cr ; t for the cho- 
sen wind temperature. The geometry of the outflow is also 
fairly unconstrained in our fits to the data from one slit; as- 
suming that the picture of [O III]-emitting clouds is correct, 
the inner half-angle of the bicone could be as small as zero, 
and the outer half-angle of the bicone could be as large as 
45°. We have chosen to r etain a bicone g eometry largely sim- 
ilar to that found in the Das et al. ( 20Q3) fits to all of their 
slits, although the outer bicone half-angle adopted here is 38°, 
somewhat larger than Das et al.'s estimate of 33° ±2°. For 
further comparison to the values in Table [2] Das et al. find a 
bicone inclination of 45° ± 5°, and an inner bicone half-angle 
of 15° ±2°. 

Using these best-fit values for the wind, we derive a mass 
outflow rate in the wind (for the case where no = 2 x 10 9 cm" 3 ) 
of M w ind = 6.6 Cf jW ind Mq yr" 1 where Cf iW i n d is the global cov- 
ering fraction of the wind. For the wind parameters found 
above, Cf, win d = 0.2, implying M wind =1.3 Mq yr" 1 . 

As mentioned previously and as indicated by the organiza- 
tion of Table 13 we can constrain only the ratios " wil " 1 and 



"wind 



A'h.ISM t 



In addition, T and n aunc h are degenerate as well; a 
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TABLE 2 

Inferred Parameters for the NLR in NGC 4151 Using Wind/Cloud Model 



Quantity 



Value 



Range Permissible for Reasonable Fit 



Wind Temperature, T 
Launch Radius, ri aunc h 
ISM Gas Radius, tism 



W H,ISM cloud 
;i wind 



^H,cloud 

Angle of Inclination of Outflow Cone 
Inner Cone Opening Angle 
Outer Cone Opening Angle 



3 X 10 6 K 
0.1 pc 
100 pc 
1.1 x 10~ 21 

6.5 x nr 21 

45° 
12° 

38° 



2.5 x 10 6 -3.5 x 10 6 
< 0.8 pc 
50-150 pc 
8.1 x 10~ 22 
n-21 



2.2 x ur 



3.2 X Ur 2u cm 



Outflow Slice in the Plane of the Sky 

— Outflow Slice out of the Plane of the Sky 
< Das et ah (8005) Data Points 




Outflow Slice in the Plane of the Sky 

— Outflow Slice out of the Plane of the Sky 
< Das et al. (8005) Data Points 



200 
r [pc] 

FIG. 6. — The best-fit velocity vs radius profiles for clouds embedded in 
a biconical, isothermal Parker wind (fit parameters are given in Table El. 
The blue lines show the velocity results for the clouds on the near side of 
the outflow, where the outflow bicone is nearly in the plane of the sky, or 
perpendicular to the line of sight. The red lines give the velocity for the 
component of the outflow bicone that is close to lying along the line of sight. 
The range of lines of each color track the velocities for intermediate outflow 
angles within the assumed bicone. The crosses give the observed velocities 
for the high-, intermediate-, and lowest-level [O III] flux emission lines from 
Das etal. I 2005). 



reasonable fit can be found at larger radii if the temperature 
drops cornmensurately. 

4. 1 . The Role of Rotation 

We do not include the effects of rotation in this model. Ro- 
tation in the Parker Wind would be important very close to the 
central source, but if angular momentum is conserved in the 
wind, V0 wind ~ 4 km s" 1 at r = 20 pc in this model. This is en- 
tirely negligible as the instrument resolution is ~ 40 km s -1 . 
The presence of rotation does slightly modify the position of 
the critical point; we have modified our ID wind models to 
include rotation (and verifi ed those m odels by reproducing, to 
within 0.1%, the results of Keppens & Goedbloed 1999) and 
found that rotation modifies the position of the Parker wind's 
critical point by only 3%; this is relatively insignificant next 
to the 23% change in the critical point caused by the adop- 
tion of the non-point-source gravitational potential. Thus, for 
these models, we have neglected the effects of rotation. 




200 
r [pc] 

FIG. 7. — As in Fig.|^]but velocities are those for the wind and not the 
clouds embedded in the wind. (Note the different vertical velocity scale on 
this plot vs. the velocity scale in Figurelol) 



4.2. The Embedded Clouds: Their Origin, Longevity, and 

Fate 

The clouds in this model are chiefly used to try to slow the 
apparent acceleration observed in [O III]: our main concern 
in this paper is to attempt to reproduce the observed kinemat- 
ics. However, it is interesting to ask about the origin of these 
clouds, their mass outflow rate, and the various effects of the 
clouds embedded in the thermal wind. 

First, a gener al overview of th e clouds in this model. It 
is suggested bv iDas et alJ J2005) that STIS can resolve the 
cloud sizes, so we consider a typical size of a cloud to be > 
lpc. For the ratio of « w ind/»cioud found in this model, M c i oul i ~ 
24.5 Ctdoud M Q yr" 1 where C f , c i ou d is of order 0.1 C t >i nd ~ 
0.02. This implies a final M c ioud ~ 0.5 M yr" 1 . 

4.2. 1 . Cloud Origins 

Most importantly, in this model, for the parameters that fit 
this data, clouds cannot be injected into the outflow at the 
launch point of the wind. The clouds must be injected down- 
stream (the injection in this wind happens at a radius of 20 pc), 
where the thermal wind's ram pressure is sufficient to over- 
ride the force of gravity felt by the clouds. If the clouds are 
injected into or exist in the model at lower radii, those clouds 
fall back to the base of the wind. 

We do not constrain the origin of these clouds any fur- 
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ther. However, it is reasonable to hypothesis that these clouds 
(given their size) are portions of large-scale clouds in the in- 
ner regions of the galaxy that are entrained in the hot thermal 
wind. One can also ask if the embedded clouds may form 
from thermal instabilities in the wind, but our Cloudy models 
of the wind are stable with the addition of adiabatic cooling 
(without such cooling, Fi g s . l3l and 1 1 41 sho w that such a wind 
at T ~ 10 6 K would be thermally unstable, but those figures 
include only heating by the central continuum and do not in- 
clude adiabatic cooling). In addition, if clouds were created 
from thermal instabilities in within the wind, they would share 
the velocity of the wind; the thermal wind itself accelerates 
too quickly to explain the slow acceleration of the [O III] 
clouds, so cloud formation by thermal instability would not 
explain the observations. 

4.2.2. Cloud Emissivity 

It is interesting to check whether the emissivity in the 
clouds is approximately equal to the observed luminosity in 
[O III]. In addition, is it possible that the clouds, which must 
be limited in mass in order to be dragged in the wind, cannot 
have enough mass to account for the \Q HI ] emiss ion? 

This is relatively easy to check. Das et alf d2005) report that 
the [O III] flux from clouds is of order 10~ 12 to 10~ 14 ergs 
cm -2 s" 1 ; in particular, near the midpoint of the clouds ac- 
celeration at 60 pc, the flux is ~ 1.6 x 10~ 14 ergs cm" 2 s" 1 . 
With NGC 4151 at a distance of ~ 13.3 Mpc, this implies a 
flux of order L[o unobserved ~ 3 x 10 38 ergs s. The clouds in 
this model have a total luminosity of L[om],theoty = 47r/'[o mjV 
where j[o nrj is the [O III] line emissivity (in units of erg 
s" 1 cm" 3 ster" 1 ) and V is the total emission volume. We 
estimate the volume as e|7rr 3 where e is the cloud volume 
filling factor, and r max is the outer radius in the wind. To 
estimate the maximum filling factor for the clouds, we set 

^^wind^windCf.wind ~ ^^wind^windCf.max.clouds wheie Cf.max,clouds IS 

the maximum possible cloud covering fraction. For the 
mass outflow rates inferred in the wind, and the cloud ve- 
locities in the wind, we find that Cf imaXiC i ou ds = 0.17; this 
is a marginally smaller covering fraction than the wind it- 
self. We then can estimate the cloud volume filling fraction, 
e < 0.02. For clouds in the best-fitting wind, we can esti- 
mate y'[oiii] ~5x 10" 21 erg s" 1 cm" 3 ster" 1 , which implies 
L[0 Hi], theory < 6 x 10 41 ergs s" 1 . It is somewhat more rea- 
sonable to assume that e ~ 6 x 10" 4 (lAlexander et aljfl999l) . 
which implies L\p nrj.theory ~4 x 10 40 ergs s" 1 ; for the sim- 
ple order-of-magnitude check employed here, this calculation 
shows that the clouds could indeed account for the observed 
emissivity, and shows that the cloud filling fraction is reason- 
ably small. 

4.2.3. Cloud Longevity 

It is also interesting to check the Kelvin-Helmholtz 
timescales against the acceleration timescales of the clouds to 
see whether these clouds might be unstable to either Kelvin- 
Helmholtz or Ra yleigh-Taylor instabilities. Both ti mescales 
are of order (e.g., Begelman, deKool & Sikora 1991) 

'/ 2 /fi> N X 1 / 2 

f KH ' 



Sv 



^wind 



cloud JV cloud 



<2.2xl0 19 s^ cloud ^ 



1/2 



^wind 

l 



l 

Jv 



\ 1 pc / Sv 
Meanwhile, the acceleration timescale is of order 



(16) 



n 1 mil -^PC 

— ~ 3.1 x 10 s — — 

V Vj 



(17) 



where R pc is the acceleration length scale in parsecs, and v 7 
is the velocity in units of 10 7 cm s" 1 . For the parameters of 
the clouds in our best-fit no = 2 x 10 9 cm" 3 model, ?kh > face 
for the clouds out to r ~ 30 pc, and ?kh is less than f a cc by 
only a factor of less than three, so only large clouds (on scales 
of 6 pc or more, which would be approximately the size of 
clouds that Das et al. seem to be resolving in NGC 4151, 
where 6 pc ~ 0.1")- Our chief concern in this model is in 
matching the observed dynamics, however. 

We can also investigate whether these clouds would be con- 
fined by the thermal pressure of the wind. Interestingly, al- 
though we did not attempt to satisfy this constraint in our 
models, the final cloud models are close to being in thermal 
pressure equilibrium with the thermally-driven wind. We do 
not model the thermal equilibrium of the clouds, but assuming 
a cloud temperature of r c i ou d ~ 10 4 K and a fiducial wind tem- 
perature of r w j na ~2x 10 6 K, clouds of radius ~ 0.2 pc would 
be in approximate thermal equilibrium. Clouds of larger size, 
such as the clouds apparently required here, would have pro- 
portionately smaller densities in our model (in order to repro- 
duce the velocity curves), and so have lower pressures for the 
assumed r c i out j. 

We note that it is possible to find a Parker wind solution 
without clouds that also fits the data, although not as well as 
the wind and cloud model; note in particular the data points 
that are not fit near (r,v) values of (20 pc, v ~ 35 km s" 1 ), 
(20 pc, 400 km s" 1 ), as well as those points above the r > 
100 pc velocity envelope of the wind. This solution is dis- 
played in Figure [8] with parameters listed in Table |3j note 
that we fit the wind-only model using the same geometry as 
used in the wind+cloud model. Since this wind-only model 
simultaneously results in a somewhat worse fit and similar 
parameters for the previous model, we consider the previous 
wind and cloud model the best representation of the data. 



Outflow Slice in the Plane of the Sky 

— Outflow Slice out of the Plane of the Sky 
< Das et al. (3005) Data Points 




FIG. 8. — The velocity vs distance profiles of the biconical outflow for the 
best fit model for the wind alone, without any embedded clouds. The data 
points are the same as those in Figures|S|andQ 



TESTS OF PHOTOIONIZATION CONSISTENCY: CAN AN 
ISOTHERMAL WIND EXIST IN NGC 4151? 
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TABLE 3 

Inferred Parameters for the NLR in NGC 4151 Using Wind Model 



Quantity Value Range Permissible for Reasonable Fit 

Wind Temperature, T 2.2 X 10 6 K 2.0 X 10 6 -2.4 X 10 6 
Launch Radius, launch 

0.1 pc <0.8pc 

ISM Gas Radius, r\sm 1 00 pc 50 - 2 1 pc 

., " wind — 1.1 x 10~ 21 cm" 1 7.4 x 10~ 22 -1.4 x 10~ 21 cm" 1 

W H,ISM cloud 

Angle of Inclination of Outflow Cone 45° 

Inner Cone Opening Angle 12° 

Outer Cone Opening Angle 38° 



Having defined an optimum temperature to fit the obser- 
vations (see Table |5Jl, we now test whether the temperatures 
required in the isothermal models can be achieved via pho- 
toionization from the central AGN. We also check for how 
the temperature varies as a function of radius in the wind, to 
test if the wind does indeed have T(r) decreasing less quickly 
than 1/r on small scales and T(r) at least consta nt on large 
scales, as required by the Parker wind model (see 32. 1> . Note 
that for the initial photoionization tests we present, the run of 
density vs. radius from the isothermal Parker- wind model is 
used to set the density at each radius in the wind, so these 
models are not completely self-consistent. We will then con- 
clude by testing an "X-ray rich" continuum to see if this can 
help explain the observations, and then finally we will test a 
non-isothermal Parker wind (with a photoionization-derived 
T(r) profile) against the observations. 

5.1. Photoionization Results 

In order to look for a physically consistent thermal wind 
model that matches the earlier isothermal model and the ob- 
servations, we ask (1) whether a roughly isothermal temper- 
ature structure can be achieved in an AGN wind, (2) whether 
the temperature is high enough to roughly match the temper- 
ature required in the isothermal wind model, and (3) whether 
T(r) decreases less quickly than M(r) / r. 

To answer the first two questions, we calculate the temper- 
ature as a function of radius for winds with a range of initial 
densities, «h,o> displayed in Figure|9] The temperature vs. ra- 
dius for each «h,o is calculated by taking the n(r) profile from 
the isothermal wind model and running Cloudy photoioniza- 
tion simulations at a range of r adii g iven the central contin- 
uum of NGC 4151 (defined in 33. 1> . This figure shows that 
only for very large «h,o is the wind approximately isother- 
mal. The only initial density in Figure |9] that leads to an ap- 
proximately isothermal wind is hh,o = 3.3 x 10 10 cm" 3 . This 
density is so high that the resultant wind temperatures rise 
only to T e max ~ 10 5 K. Such a temperature is far too low to 
explain the observed velocities in NGC 4151 which required 
T e = 3 x 10 6 K. Lower densities («h,o 3.3 x 10 9 ) are also 
not possible in this model as T(r) always drops significantly 
at large distances. This is simply due to adiabatic cooling. 
With adiabatic cooling switched off, the temperature in the 
wind very quickly achieves T e w 3 x 10 7 K and remains at that 
constant temperature over the rest of the outflow. Even with 
adiabatic cooling included, however, Figure|9]also shows that 
«h,o ~2x 10 9 cm" 3 can achieve approximately the required 
temperatures over at least a limited range of radii in the out- 
flow; the drop in T(r) at larger distances shows that it may 
still be problematic, however. 



Next, we focus on how quickly the temperature can change 
with radius and still allow a Parker wind to form. Figure 
displays the exponent of a piecewise power law that fits the 
temperature profile shown in Figure |9] Again, in order for 
a physically consistent Parker wind to exist in the core of 
NGC 4151, T(r) must drop less slowly than 1/r for small 
distances (r < 2 pc), and then, where M(r) increases almost 
linearly with radius, T(r) must be, at least, approximately 
constant: it must have a power law index > 0. This is true 
for the high density run with «h,o = 3.3 x 10 10 cm" 3 , but as 
already shown in Figure [9] that wind does not have tempera- 
tures high enough to yield the observed velocities. The lower 
density trials meanwhile show extraordinarily large changes 
in r e with radius. For n < 3.3 x 10 8 cm" 3 , T(r) drops be- 
low a power law index value of zero at r < 0.7 pc. In con- 
trast, acceleration of the isothermal Parker wind must con- 
tinue out to a few times the Parker wind's critical point dis- 
tance, which is at r cr i t ~ 1 pc for the isothermal wind model 
shown in Figure [6] (non-isothermal wind models have critical 
points even further out, since much of the wind's acceleration 
occurs with T(r) < 7i sot hermai)- However, the temperature is 
already dropping precipitously at radii even smaller than r cl -i t 
for «h,o < 3.3 x 10 8 cm" 3 , so that range of densities would not 
yield a Parker wind. 

The test case of «h,o = 2 x 10 9 cm" 3 is closer to having the 
required temperature (at least at its maximum temperature), 
and is closer to the desired T(r) profile, and so more careful 
consideration is justified in this case. As shown in Figure [9] 
this model does, if only in a limited range of radii, achieve 
approximately the temperature of the best-fitting isothermal 
model. Figure [K)] shows that the temperature does not drop 
appreciably until r ~ 2 pc. We would expect that since accel- 
eration in the wind is required out to a few parsecs, this tem- 
perature profile would also not lead to a viable Parker wind 
model. But, as the temperature profile is close to satisfying 
the requirements on T(r), we build a numerical Parker wind 
model for generic T(r) from our photoionization simulations 
that finds, self-consistently, the critical point for that particu- 
lar temperature law. Applying this model to T(r) for the case 
of «h,o = 2 x 10 9 cm" 3 , we found that the wind decelerates 
beyond r = 8 pc, and so would not fit the required v(r) pro- 
file. Testing to see if this «h.o = 2 x 10 9 cm" 3 profile results 
in clouds velocities that might fit the data, we apply the same 
geometry and dragparameters as used in the previous best-fit 
model (see Table and present the model vs. the data from 
iDas et alJ J2005I) in Figure The model clearly cannot fit 
the observed velocities, or their variation with distance. 

If we instead use the T(r) profile derived from the «h,o = 
3.3 x 10 8 cm" 3 wind model, we again find a wind that deceler- 
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ates too quickly (here, the wind starts decele rating at r = 4 pc) . 
Attempting to fit this model to the data from Das et al. (2005) 
results in a much better fit (see Figure [121 than in Figure ITTI 
but still misses the high velocity data points shown in the fig- 
ure. Modifying the drag term in our wind equation does not 
improve the fit: decreasing the drag term to fit the higher- 
velocity points simply increases the outflow velocities, there- 
fore doing an increasingly worse job of fitting the data points 
near (r, v) = (70 pc,350 km s -1 ). Also, modifying the geom- 
etry of the outflow has no significant impact on the fit: the 
half-angle of the bicone is already 38° here; increasing that to 
a maximum value of 45° helps to include some of the higher 
velocity data points, but cannot make up for the discrepancy 
relative to the isothermal model. Still, this model represents 
the closest fit to the data achievable with a more realistic T(r) 
profile; there is no quantitative way to distinguish this fit from 
the previous isothermal best-fit, but it is perhap s enc ouraging; 
this model will be examined in more detail in 35.21 Unfortu- 
nately, all of the rest of the generic T(r) models cannot come 
close to reproducing the data (as can be seen, for example, in 
Fig.O- 
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FIG. 9. — Tc in the wind as a function of radius in the Parker wind. As 
shown in the figure, both the rise in T c as a function of radius and the max- 
imum T e achieved are strongly dependent on the initial density. Only winds 
with «h,o ^2x 10 9 cm" 3 achieve temperatures of T ~ 3 X 10 6 K required by 
the isothermal Parker wind model. Note also the slower increase of T t with 
radius as the density increases. 



5.2. How does a self-consistent T(r) affect the Parker wind 
simulations? 

Up until this point, we have used the n(r) profile from an 
isothermal Parker wind model as input to a set of photoion- 
ization models, only rescaling the base density to test different 
initial densities. It is important to check if the fitting results 
improve with a more self-consistent T(r), n(r), and v(r) pro- 
file instead of using the isothermally-derived n(r) profile to 
calculate T(r), as before. 

We test this possibility by examining the T(r) profile from 
the no = 3.3 x 10 8 cm 3 simulation (shown in Fig.|9}. This is 
then input into our non-isothermal Parker wind model, and 
the resulting acceleration profile and n(r) profile is computed. 
This n(r) profile is then used as the basis for another itera- 
tion of Cloudy simulations, to compute the resultant T(r). We 
ask whether the result of this iteration on the T(r) structure 
will result in a temperature profile more suitable for a large- 
scale Parker wind. However, the resulting T(r) profile only 
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FIG. 10. — A piecewise power-law fit to the temperature data in Figure l9l 
with the power law, a defined via T oc r a for radii between the photoion- 
ization simulation results. Physically consistent Parker winds in the NLR of 
NGC 4151 would require T(r) dropping less quickly than \jr out to r ~ 2 pc, 
and T(r) at least constant out to r ~ 10 pc in order for the wind to acceler- 
ate. As can be seen here, this is very difficult for any of the winds (given the 
dominance of adiabatic cooling at large radii) to satisfy. 
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FIG. 11. — As in Figure [5] but with a Cloudy-derived T{r) profile for 
«H,0 = 2 X 10 9 cm -3 . The wind, with a maximum velocity of only ^365 
(the cloud velocities, shown here, do not even approach that velocity), 
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and decelerating for r > 8 pc, cannot drag the clouds to anywhere near the 
observed velocities. 



results in a worse fit to the data: the rise in the initial T(r) 
profile at r < 1 pc causes the wind model to accelerate from 
the disk faster than the fiducial isothermal wind model. This 
increase in acceleration yields a wind that must be launched 
from larger radii (r ~ 0.2 pc instead of 0.1 pc) and also has 
a faster dropoff in density with height, which then yields a 
T(r) curve closer to the lower-density curves in Figure |9] 
the self-consistent T(r) only results in a less isothermal wind 
model whose temperature decreases more quickly with radius 
beyond r « 1 pc. Therefore, attempting to include a self- 
consistent temperature profile in the Parker wind model only 
results in a more rapidly increasing (and then more rapidly 
decreasing) T(r) profile, which results in a wind model that 
has a T(r) profile even less s uited for accelera ting a wind that 
might fit the observations o f lDas et al.ld2005l) . So, in the end, 
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FIG. 12. — As in Figure|6] but with a Cloudy-derived T(r) profile for 
n H .o = 3.3 X 10 s crrT 3 . 



including such a T(r) profile in the Parker wind also cannot 
explain the observed acceleration in NGC 4151. 

5.3. Does an X-ray "rich" continuum help? 

To test the limits of our model of the central continuum, 
we now adopt a relatively X-ray-bright continuum to check 
if such a continuum would heat the gas to higher temper- 
atures and modify our results. Displayed in Figure 
this continuum was primarily modified in the high-energy 
X-ray regime to m ore closely follow the observations of 
Edelson et alJ i 19 961) and the inferred high-energy spectrum 
of Alexander et al. (1999); we continue to include the UV 
absorption, as before. This increase in X-ray brightness in- 
creases Compton heating within the gas. Meanwhile, the 
low-energy continuum in the SED has been significantly cut 
off below 1 eV in order to reduce Compton cooling by low- 
energy ph otons as much as possible. This was done, for in- 
stance, in Rrolik ( 1999) to account for the possibility that 
lower-energy photons are emitted at much larger scales and 
so might not act to cool the outflow in the immediate vicinity 
of the central black hole. (In reality, for the large-scale winds 
tested here, the question of the origin of the low-energy pho- 
tons is less important since they almost certainly play a role 
on scales of 10 pc. However, to test the limits of the input 
continuum model, we attempt such a continuum anyway.) 

Adopting this continuum does not greatly modify the 
Compton temperature of the gas. The radiative equilibrium 
curve for this continuum is shown in Figure 1141 for this 
continuum, r Co mpton = 4.2 x 10 7 K as opposed to r Co mpton = 
3 x 10 7 K for the original continuum. 

As one might infer from that relatively small change in the 
Compton temperature of the gas, adopting this continuum for 
calculating T(r) does not result in significant changes to our 
results. We also calculated a self-consistent Parker wind after 
inputing this continuum into our photoionization simulations, 
which yielded no appreciable difference from the tests with 
the original SED. Overall, adopting this SED does not signif- 
icantly alter any of our conclusions. An X-ray "rich" contin- 
uum for the central SED of NGC 4151 cannot help a Parker 
wind model to explain the observed velocity trends. 

6. CONCLUSIONS 
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FIG. 13. — An alternate X-ray "rich" incident continuum for NGC 4151, 
plotted with the same data as in Fig.[2| For this continuum, w e explicitly set 
the high-energy X-ray continuum to fit Edelson et al. 1 1996), as opposed to 
the previously adopted continuum, which adopted the Rraemer et al. 1 2000) 
power law in the X-rays. This led to slightly higher X-ray flux. At the same 
time, the continuu m below 1 eV was suppressed to lower Compton cooling 
in the outflow (as in lKrolrkll9^9T) . 
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FI G. 14. — The photoionization equilibrium curve for the SED shown in 
Fig. 1131 As can be seen by comparing to Fig. [3] 7c om p ton does not change 
appreciably. 



We have developed models of Parker winds accounting for 
the effects of the radially varying M(r), adiabatic cooling, and 
drag effects between the wind and embedded clouds as well 
as drag forces between the wind and an external medium. We 
then tested these models against the conditions that exist in 
the center of one particular local AGN, NGC 4151. We have 
shown that simple thermal winds cannot explain the slow ac- 
celeration that has been inferred for the NLR of NGC 4151. 
The strongest constraint against self-consistent Parker winds 
is that adiabatic cooling on the large scales of these outflows 
leads to temperatures in the wind that decreases strongly as a 
function of r. With such a temperature profile, Parker winds 
cannot be launched and accelerated to the observed velocities 
at the observed distances. Varying the initial density in the 
wind either yields temperatures that drop quickly with dis- 
tance (for «h,o ^ 10 8 cm -3 ) or temperatures that are too low 
to reproduce the observed velocities (n H ,o > 10 9 cm -3 ). The 
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intermediate case of «h,o ~ 3 x 10 8 cm -3 comes closest to 
reproducing the data, yielding a temperature profile that we 
have checked by developing a numerical Parker wind model 
for general T(r). For the temperature profile shown in Fig- 
ure |9l the wind starts decelerating already at r ~ 4 pc; a ra- 
dius much smaller than the 100 pc out to which acceleration 
is inferred in NGC 4151. We find that plo tting the re sultant 
predicted v(r) against the data from D as et alJ (l2005) shows 
that this wind model cannot reproduce the range of [O III] 
velocity points (see Fig. I12> as well as the isothermal wind 
model (which cannot be obtained for NGC 4151). Further, at- 
tempting to develop a self-consistent T(r) profile for this wind 
results in a temperature profile even further from isothermal 
which correspondingly fails to fits the data. In short, we find 
no self-consistent T(r) profile for a large-scale Parker wind 
that can explain the observed ranges of velocities. 

If thermal winds are not driving the observed v(r) profile, 
how can it be explained? Radiative driving is most likely not 
the source of the acceleration: we have tested simple mod- 
els of radiative acceleration on dust, and such models achieve 
terminal velocities at distances of order the launch radius - 
much closer to the launch point than the observed velocity 
profiles. Magnetic winds would also achieve their terminal 
velocities on similar scales (although perhaps clouds dragged 
in the magnetic wind might yield the slow acceleration seen). 

The difficulty in fitting this data with almost any known 
wind model leads us to question whether the observed veloc- 
ity law is actually a wind accelerating with radius. Perhaps the 
observed kinematic profile is the result of an interaction of a 
wind with the surrounding medium. An interesting alternative 
interpretation is suggested in the work of M atzner & McKeel 
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(1999), which showed that rather generic hydromagnetic 
winds, expanding into gas surrounding a young stellar object 
at the wind's terminal velocity, would yield v oc r, or a "Hub- 
ble law"-type flow. This is a promising theory for the acceler- 
ating portion of the outflow, but it seems unclear how the ob- 
served slow deceleration of the flow in NGC 4151 could then 
be explained in the same picture. Still, it is important to point 
out that given the inability of a wide range of wind models to 
explain the observed v(r) profile, perhaps these observations 
are not revealing acceleration of a wind, but instead the in- 
teraction of an already accelerated wind with the surrounding 
medium. 
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